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ABSTRACT 

We show how, based on considerations on the observed form of the galaxy 
two-point spatial correlation function ^(r), a very simplified - yet surprisingly 
effective - model for the linear density fluctuations power spectrum can be 
constructed. We flrst relate the observed large-scale shape of ^(r) to a power- 
law form for the power spectrum, P{k) oc /c~^'^. For a plausible value of the 
bias parameter b = l/cs — 1-8, one has {dp/p)rms ~ 1 at r ~ 3.5 h~^ Mpc, 
suggesting that the change of slope observed in ^(r) around this scale marks 
the transition between the linear and nonlinear gravitational regimes. Under 
this working hypothesis, we use a simple analytical form to flt the large-scale 
correlations constraints together with the COBE CMB anisotropy measure- 
ment, thus constructing a simple phenomenological model for the linear power 
spectrum. Despite its simplicity, the model flts remarkably well directly esti- 
mated power spectra from different optical galaxy samples, and when evolved 
through an N-body simulation it provides a good match to the observed galaxy 
correlations. One of the most interesting features of the model is the small- 
scale one-dimensional velocity dispersion produced: aid = 450 Km s~^ at 0.5 
h~^ Mpc and aid = 350 Km s~^ for separations > 2 h~^ Mpc. 



1. INTRODUCTION 

Cosmological models for the formation of the large-scale structure of the 
Universe require as a key ingredient a speciflc form for the power spectrum of 
density fluctuations P{k) at recombination. This is usually specifled in terms 
of a primordial spectrum together with a transfer function which describes the 
subsequent evolution of fluctuations depending on the physical scenario. A 
paradigmatic example is provided by the Cold Dark Matter (CDM) model (e.g. 
Blumenthal et al. 1984; White et al. 1987). 

Conversely, one might ask, independently of any a priori physical model, 
whether a speciflc shape for the power spectrum is implied by the available 
data, i.e. adopt a pLeiiomeuoiogicai approach (e.g. Kashlinsky 1992; Scaramella 
1992; Taylor & Rowan-Robinson 1992). Indeed, observations have reached a 
stage in which it is possible to constrain directly the power spectrum on linear 
scales, where {dp/p)rms < 1 and the recombination shape of P{k) should be 
preserved. On one side galaxy redshift surveys have allowed reliable estimates 
of the two-point correlation function ^(r), the Fourier counterpart of P{k), up 
to separations ~ 30 h~^ Mpc (de Lapparent, Geller & Huchra 1988; Guzzo et 
al. 1991, hereafter G91; Loveday et al. 1992), and direct estimates of P(A;) itself 
(Baumgart & Fry 1991; Peacock & Nicholson 1991; Vogeley et al. 1992, here- 
after V92; Jing & Valdarnini 1993, hereafter JV93; Fisher et al. 1993, hereafter 
FDSYH). On much larger scales, the discovery by the COBE-DMR experiment 
of temperature fluctuations in the microwave background radiation (Smoot et 
al. 1992) has produced a further fundamental constraint on the very long- 
wavelength amplitude of P{k). Among the outcomes of this relative wealth 
of results is the fact that the standard flat CDM model (although successful 
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on several grounds), cannot consistently reproduce all the observations. The 
difficulties arise essentially from a too-high ratio of small-scale to large-scale 
power. Several modifications to the standard CDM scenario have been explored 
in the attempt to overcome these problems . Among these, open models with 
non-vanishing cosmological constant (Efstathiou et al. 1990), models with tilted 
primordial index (Adams et al. 1992, hereafter A92; Cen et al. 1992), or mixed 
(CDM plus HDM) models (Davis, Summers & Schlegel 1992 and references 
therein). Rather than entering the debate on which physical model is closer to 
the observed behaviour, in this Letter we prefer to keep to the pure phenomeno- 
logical approach to reconstruct a "minimal" model for the linear P{k). As we 
shall see, when we close the circle, despite its naivety the model proves to be 
remarkably effective in accounting for a number of observations. 

2. A PHENOMENOLOGICAL POWER SPECTRUM 

The most natural way to construct a phenomenological model for the linear 
P{k) would be to start from its most recent direct determinations (as e.g. V92 
or JV93), try to 'clean' the observed shape from the expected non-linear ef- 
fects, and then add further constraints from other observables. Here, however, 
we choose to start in an apparently less straightforward way, i.e. from the ob- 
served two-point correlation function ^(r). The reason is that we first want to 
show how the large-scale drop-off of ^(r) contains important information and, 
although not a power-law, can be related to a specific, simple shape of P{k) 
on intermediate scales. We stress that the phenomenological model obtained 
from this heuristic procedure will be fully justified only a posteriori, when it will 
possibly satisfy the main observational constraints. 

Let us therefore consider the observed spatial correlation function estimated 
from redshift survey samples. This is different from the true (real space) one, be- 
cause on small scales large velocity dispersions within clusters ('Fingers of God') 
depress its amplitude, while on large scales coherent motions amplify it accord- 
ing to linear theory prescriptions (Kaiser 1987). These effects can be properly 
accounted for by estimating ^(rp,7r), (e.g. Fisher et al. 1993b), i.e. decompos- 
ing the pair separation vector along the directions parallel and perpendicular to 
the line of sight. A practical, approximate correction of the small-scale distor- 
tions can however be performed by collapsing the cluster 'fingers' into a region 
corresponding to their statistically expected spatial size. Using this method, 
G91 obtained from the Perseus-Pisces (PP herafter) redshift survey, a corre- 
lation function reproducing the canonical real-space shape for optical galaxies 
oc r~^'^ (Davis & Peebles 1983), shown in fig. 1. Although this function is still 
in redshift space at large separations, we prefer to use the notation ^(r) instead 
of the conventional ^(s) to underline the correction to real space operated on 
small scales. In the same figure we plot (open circles) the more recent estimate 
(in pure redshift space) from the extension of the CfA redshift survey (V92). 
Note the the small-scale depression of the latter correlation function and the re- 
markable agreement of the two samples for r > 2 h~^ Mpc. In particular, both 
samples are very well described for separations larger than r^ ~ 3.5 h~^ Mpc 
by the simple law ^(r) = [(r/20)~'^'^ — 1] (dashed line), where r is in h~^ Mpc. 
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The consistency of the two surveys further demonstrates that the original find- 
ing of G91 cannot be ascribed to pecuharities in the Perseus-Pisces region, and 
that Tb ~ 3.5 h~^ Mpc seems to represent a physicaUy meaningful scale (see 
also Calzetti, Giavalisco & Meiksin 1992). 

By Fourier transforming the above expression for ^(r) - assuming negligible 
correlations for s > 20 h~^ Mpc - we can have a first rough indication of the 
shape of the corresponding power spectrum: P{k) oc /c~^'^. Although ^(r) is 
not rigorously a power law at large separations, the corresponding P{k) still 
seems to be. More accurately, if we adopt a simple functional form as used by 
Peacock (1991) to fit the APM angular function, 

P(k) = f ^'° „-,. (1) 

i + (4 

then we can match the observed correlations for s > ri, with a ~ —2.2, 
kc — 0.08 h Mpc~^ and virtually any value of n in the range [0, 1] (which 
has very little effect on ^ for A < 27: /kc)- At this point we have essentially 
two ways to constrain the primordial index n. The first is simply to assume an 
Harrison-Zel'dovich primordial spectrum, i.e. n = 1. In this case, considering 
the microwave background rms temperature fiuctuation on 10° scales measured 
by COBE, aT(10°) = [1.085 ± 0.183] x 10"^ (Smoot et al. 1992), we obtain a 
normalization corresponding to a bias factor b = l/cxs — 1-4, where ag is the 
mass variance in a top-hat window of 8 h~^ Mpc radius. For clarity, note that 
we are using here the conventional definition of 6 = as{gal) / as{mass), together 
with the Davis & Peebles (1983) result that gives as{gal) ~ 1. One may define 
a bias factor for bright galaxies, 6*, compatible with the correlation functions of 
fig. 1, considering that, once the redshift space amplification is corrected using 
Kaiser (1987) formula, is a*s{gcil) — 1.3 and therefore b* ^ 1.36. 

A completely different route to get a plausible value of the (n, b) pair is to 
note that if 6 = 1.8 — 2, a range which turns out to be justified in several 
respects for optical galaxies, then we have a^; ~ 1 for i? ~ 3.5 h~^ Mpc = r^. 
In other words, in such a case the change of slope in the observed correlation 
function is located right around the scale where {Sp/p)rms ~ 1, i-e. coincides 
with the expected transition from the linear to the nonlinear clustering regimes. 
This has the important consequence that larger scales should still be evolving 
in a quasi-linear fashion and could therefore constrain directly the shape of 
the initial spectrum. Using this kind of bias value, and considering again the 
COBE constraint at large wavelengths we obtain n ~ 0.75 if we neglect the 
gravitational wave contribution to the Sachs- Wolfe effect (Lucchin, Matarrese & 
MoUerach 1992; A92). We note that in this case all the important parameters of 
the spectrum, i.e. b ~ 1.9, a ~ —2.2 and n ~ 0.75 are globally consistent with 
those of a correspondingly tilted CDM spectrum (Cen et al. 1992; A92). This is 
evident in fig. 2a, where we compare our euristic fit with CDM models of varying 
primordial index (n = 0.7, 0.8, 0.9, 1), adopting the Bardeen et al. (1986, BBKS) 
CDM transfer function. It is implicit in this comparison that we are considering 
our large-scales fit as representative of the linear power spectrum, in particular 



that any residual redshift space effects are limited to the Kaiser amplification. 
This is quite a strong assumption and might be questionable in several ways 
(e.g. FDSYH; Bahcall, Cen & Gramann 1993). However, let us assume here 
as a working hypothesis that expression (1), with a = —2.2, n = 0.75 — 1, 
b = 1.4 — 1.9, represents a good linear phenomenological description of the 
true P{k). We shall call it the Phenomenologically Induced Model (PIM). The 
remainder of the paper will concentrate on comparing it with the observations, 
and possibly narrowing the constraints on its normalization and primordial 
index. 

3. LINEAR AND NONLINEAR TESTS 

The first test is obviously to check whether our phenomenological form is a 
good representation of the 'general' behaviour of clustering, as quantified by 
power spectra measurements on independent samples. To this end, in fig. 2b 
the PIM, with n = 0.75 and n = 1 and normalized to fit the observed large- 
scale ^(r) of the PP sample, is compared to the direct estimate of P{k) from the 
CfA survey (V92). This comparison is particularly significant, since both sur- 
veys are constructed starting from the same photometric material (the Zwicky 
catalogue), and both the PP ^(r) and the CfA P{k) are estimated on volume- 
limited samples with M brighter than ~ M*. For these reasons we prefer to 
avoid including in the figure measures of P{k) for, e.g., radio galaxies (Peacock 
& Nicholson 1991), or IRAS galaxies (FDSYH; JV93), which inevitably require 
an arbitrary amplitude renormalization to be compared with the optical data 
(see fig. 1 of V91 for an indirect comparison). The agreement of the PIM with 
the CfA P{k) is evidently very good. Note that there is no renormalization 
between the curve and the data points. 

What seems to be worrying in fig. 2b is that the data points continue to 
follow the PIM curve - which is supposed to be describing the linear P{k) - 
also for for k > 0.3 h Mpc~^, where nonlinear effects certainly start to become 
significant. To understand this, we have Fourier-transformed the PP ^(r) of 
fig. 1, which is corrected for the main nonlinear distortions produced by the 
'Fingers of God'. The result is plotted as a dashed line in in fig. 2b, showing 
that, as one expects from the small-scale r~^'^ shape of the correlation function, 
P{k) rises to ~ /c~^'^ as consequence of nonlinear clustering. It is interesting 
to notice, therefore, that there seems to be a kind of 'conspiracy': nonlinear 
clustering enhance correlations at small-scales, but when these are viewed in 
redshift space, the global effect of small-scale velocity dispersions seems to bring 
the spectrum back to a shape similar to the purely linear one. This is the same 
effect visible in the correlation function (fig. 1): when the 'finger of God' effect is 
not corrected (CfA), ^(r) is well approximated by a single law oc [(r/20)~'^'^ — 1] 
from small to large scales. 

The other important test is to check how the PIM behaves when nonlinear 
gravitational evolution sets in. To this end, we have performed an N-body sim- 
ulation, concentrating our attention on two of the most crucial observational 
tests: the two-point angular correlation function w{9), and the one-dimensional 
galaxy velocity dispersion aid- We do not perform here a detailed investigation 
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of redshift-space clustering (as e.g. in FDSYH), which is left for a more com- 
prehensive paper involving further simulations with extended dynamic range. 

For our specific aim, this first simulation was planned as a compromise be- 
tween the necessity of a large box size, required by the considerable amount of 
large-scale power present in the PIM spectrum {Lhox > '^'^/ktum ~ 130 h~^ 
Mpc), and the need of a good spatial resolution for not to underestimate aid- 
We therefore used a P"^M code to integrate the motion of 64"^ particles within a 
180 h~^ Mpc sided box, with 160^ grid points and a resolution of 0.3 h~^ Mpc. 
We used Oq = 1 and Ho = 50 km s~^ Mpc~^. Initial conditions were described 
by a Gaussian random field, and the 'galaxies' were selected using the BBKS 
prescription, as in White et al. (1987). 

A further important reason for performing an N-body test is to possibly put 
some more stringent constraints on n and h. To this end, we adopted the 
following procedure. We did not consider the COBE normalization to fix the 
present epoch, but leave clustering evolve until the small-scale shape of ^(r) for 
'bright galaxy peaks' (defined with a threshold corresponding to I'sim = 1.36 ^ 
2.3) in the simulation matched the observed one. At this point h = l/crg ^ 
1.8, and the corresponding n required to match COBE with the PIM form is 
n ~ 0.85. Fig. 1 shows the two-point correlation function for the simulated 
galaxies (solid line), calculated in real space as in White et al 1987. In order 
to be compared with the PP estimate, ^(r) has been multiplied by the proper 
linear amplification factor (Kaiser 1987), on the scales where {6p/p)rms < 1- 
The match with the observations is remarkable in both shape and amplitude 
over a wide range of scales, a result which seems to support nicely the idea 
that the change of slope at r^ ~ 3.5 h~^ Mpc tags the transition scale to full 
nonlinearity. 

A crucial test for the model is provided by the angular correlation function 
w{9), which is by definition independent from any redshift space effect. We 
estimated w{9) by projecting ^(r) of the simulation through the Limber equa- 
tion, adopting a selection function and normalization corresponding to the Lick 
depth (see Peacock 1991). Since at large separations (r > 15 h~^ Mpc) the 
single simulation performed shows slightly less power than the linear model, on 
these scales we used ^(r) from the latter one, smoothly joined to the numerical 
result at smaller scales. The final result is shown in fig. 3, compared with the 
APM (Maddox et al. 1990), and EDSGC (Nichol & CoUins 1993) data. The 
PIM spectral shape does still seem to lack some power at very large angular 
separations. This might either suggest that the spatial estimates of ^(r) and 
P{k), on which the PIM is based, still suffer from the limited sample size, or 
be related to possible sistematics in w{9) at these large angles. 

The other main aim of the simulation was to check the small-scale velocity 
dispersion produced by the PIM. We computed aid{r) for the dark particles as 
o'ldi'f) — l/\/3(|vi — V2|)^/^, where vi and V2 are the peculiar velocities of the 
two particles with separation r. We obtained a"id(0.5 h~^ Mpc) ~ 450 Km s~^, 
rapidly decreasing to ~ 350 for r > 2 h~^ Mpc. This compares interestingly 
with the observed aid — 320 Km s~^ on 1 h~^ Mpc scales (Fisher et al. 1993b). 

In summary, we have shown that an extremely simplified toy model with a 
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few specific features, in particular a steep slope ~ A;~^'^, a turnover around k ~ 
0.05 h Mpc~^ and a biasing 6 ~ 1.8, is capable of reproducing remarkably well 
some of the main properties of galaxy clustering. This supports the soundness 
of the basic assumptions made for its construction, in particular the idea that 
fiuctuations on scales larger than ~ 5 /i~^ Mpc are still weakly affected by 
nonlinear evolution. Clearly, more extended tests are required, in particular 
direct comparison of real and redshift space properties from several simulated 
catalogues. Also, we have not touched here the issue of which physical scenarios 
might be able to produce a spectrum at recombination with this form. However, 
it seems that this first phenomenological investigation has at least underlined 
a few features which seem to be unavoidable for any physical model aimed at 
explaining the observed large-scale structure. During the revision of this paper, 
we became aware of a work by Baugh & Efstathiou (1993), who accurately 
determine P{k) in true real space by de-projecting the APM w{9). Their result 
fully agrees with our claim that the power spectrum has a steep range ~ k~'^ 
for 0.08 < k < 0.3 h Mpc~^ which is not the result of redshift-space effects, 
and raises to /c~^'^ for k > 0.3 h Mpc~^ due to the onset of nonlinear clustering 
(cf. our fig. 2b with, e.g., their fig. 12). 
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Figure captions 

Figure 1. 

Two-point correlation functions from the PP (G91) and CfA (V92) redshift 
surveys. The CfA estimate is in pure redshift space, while the PP sample has 
been statistically corrected for the small-scale 'Finger of God' effect. Error bars 
are 1-a bootstrap errors. The dashed line shows the curve ^(r) = (r/20)~*^'^ — 1, 
while the solid line is the result of the P'^M N-body integration of the PIM. 

Figure 2. 

a) The power spectrum of the phenomenologically induced model (PIM), with 
primordial index fixed to n = 0.75, normalized to COBE and compared to flat 
CDM models with varying primordial index (n = 0.7, 0.8, 0.9, 1). 

b) The PIM power spectrum normalized to match the PP galaxy correlations 
of Fig. 1 (solid line) for the two cases n = 0.75 and n = 1, compared to the 
direct estimate of P{k) from the CfA survey by V92 (circles). The dashed line 
shows the result of Fourier transforming the observed PP small-scale ^(r). 

Figure 3. 

Angular correlation function w{9) calculated from the N-body integration 
of the PIM (solid line), compared with the APM (dots) and EDSGC (circles) 
data. The EDSGC la error bars show the large uncertainty existing in w{9) 
for ^ > 10°. 
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